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I.  INTRODUCTION 


A.  The  Discovery  of  the  Zodiacal  Dust  Bands 

One  of  the  major  discoveries  of  the  Infrared  Astronomical  Satellite  (IRAS)  was  of  three 
parallel  bands  of  dust  roughly  straddling  the  plane  of  the  ecliptic  (Low  et  al.  1984).  Zodiacal 
dust  bands  were  first  noticed  in  the  IRAS  data  as  a  pair  of  symmetrically  placed  bumps  of 
12  and  25  micron  emission  superimposed  on  the  smooth  zodiacal  background  (Figure  1). 
The  same  emission  could  be  seen  in  the  60  and  100  micron  data  at  a  much  lower  intensity, 
indicating  a  fairly  high  temperature  (about  200K)  for  the  emitting  material. 

During  its  All-Sky  Survey  the  IRAS  telescope  scanned  the  sky  in  circles  of  con¬ 
stant  solar  elongation,  its  nearly  polar  orbit  precessing  by  ~  1°  a  day  in  order  to  remain 
above  the  terminator  on  the  earth.  The  entire  sky  can  be  mapped  over  a  period  of  six 
months  in  this  way.  More  specifically,  for  the  first  two-tiiirds  of  its  ten  month  life,  IRAS 
used  elongations  between  80°  and  100°  to  scan  nearly  the  entire  sky  4  times  (Neugebauer  et 
a/.1984;  IRAS  Explanatory  Supp.  1988).  Each  scan  had  a  width  of  0t5,  and  was  shifted  in- 
longitude  by  about  0!25  on  the  subsequent  orbit,  observing  by  overlapping  scans  the  same 
location  twice  after  ~  103  minutes.  This  was  called  an  “hours-confirmed”  observation  or 
HCON.  The  initial  observing  strategy  was  to  allow  the  telescope  to  map  a  section  of  sky 
in  this  fashion  by  slowly  changing  its  solar  elongation  for  about  a  week  (HCON  l),  after 
which  the  telescope  was  repositioned  and  the  same  section  of  sky  was  observed  a  second-  . - 
time  (HCON  2). 

As  the  sky  survey  unfolded,  the  emission  bumps  were  found  to  be  distributed 
more  or  less  completely  around  the  solar  system.  The  initial  impression  was  that  two  disks 
or  belts  of  material  had  been  found,  presumably  in  the  inner  solar  system  because  of  their 
geometry  and  high  color  temperature.  Such  independent  bands  would,  at  some  longitude, 
cross  each  other  and  the  ecliptic.  However,  as  the  mapping  became  more  complete  the 
emission  was  found  to  form  parallel  bands  above  and  below  the  ecliptic  plane.  They  did  not 
cross.  This  configuration  was  difficult  to  explain  at  first,  since  the  upper  and  lower  bands 
were  thought  a  priori  to  be  composed  of  separate  groups  of  material  which  would  have  to 
cross  the  mid-plane  given  the  Sun’s  central  force.  IRAS  Science  Team  member  D.  Beintema 
suggested  the  correct  explanation:  A  band  pair  can  arise  from  a  single  distribution  of 
material  with  an  ensemble  of  orbits  which  shared  a  common  inclination,  but  whose  nodes 
were  uniformly  distributed  over  all  ecliptic  longitudes.  This  distribution  produces  a  band 
pair  because  each  individual  particle  spends  most  of  its  time  at  its  extreme  separation  from 
the  ecliptic,  travelling  roughly  parallel  to  the  plane,  much  like  a  pendulum  is  preferentially 
found  at  its  maximum  amplitude.  A  particle  spends  little  time  near  the  ecliptic  plane 
because  there  it  has  a  large  velocity  component  normal  to  the  plane. 

A  high-pass  spatial  filter  can  be  employed  to  remove  the  smooth  zodiacal  back¬ 
ground  from  the  IRAS  survey  scans,  revealing  many  details  of  the  zodiacal  dust  bands. 
This  process  was  used  to  produce  the  emission  maps  shown  in  Figure  2.  The  originally 
recognized  band  pair  (7)  is  seen  ~  10°  above  and  below  the  ecliptic.  Part  of  the  central 
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band  is  due  to  the  peak  in  the  smooth  zodiacal  emission  near  the  plane  leaking  through 
the  high-pass  filter,  but  there  is  also  clearly  another  pair  of  bands  similar  to  the  10°  bands, 
but  spaced  only  one  or  two  degrees  from  the  ecliptic  plane.  These  inner  bands  were  later 
separated  further  into  two  pairs  of  bands  (Sykes  1986). 

Further  analysis  of  the  IRAS  data  has  suggested  that  four  additional  pairs  of 
bands  may  also  exist  (Table  I)  (Sykes  1988),  and  that  these  bands  extend  over  40°  of 
ecliptic  latitude.  In  the  IRAS  skyflux  maps  the  central  bands  (a  and  0)  have  apparent 
widths  of  <  0!5.  The  other  bands,  particularly  the  7  bands,  are  several  degrees  in  apparent 
width  (Sykes  1988).  As  will  be  seen  later,  the  broad  morphology  of  most  dust  bands  has 
a  significant  impact  on  the  number  of  bands  which  are  possible  to  observe.  So,  from  the 
original  three  bands  reported  by  Low  et  al. (1984),  as  many  as  14  bands  (7  pairs)  have  now 
been  detected. 

B.  Initial  Analysis 

During  its  last  months  of  operation,  IRAS  surveyed  the  sky  a  last  time  (HCOX 
3),  but  changed  its  scanning  pattern  to  begin  at  extreme  solar  elongations  of  60°  and  120° 
and  smoothly  approach  90°  elongation  from  both  sides  and  then  reverse  the  procedure  until 
the  whole  sky  was  covered  (Figure  2b).  This  map  includes  only  72%  of  the  sky  because  the 
survey  was  prematurely  terminated  by  liquid  helium  exhaustion  aboard  the  satellite.  The 
7  bands,  in  particular,  show  clear  variation  in  their  separation  corresponding  to  changes  in 
viewing  geometry,  and  can  be  readily  understood  if  the  density  enhancements  giving  rise 
to  the  bands  are  confined  to  a  small  range  in  distance  from  the  sun,  maintain  a  constant 
linear  separation  around  their  circumference,  and  lie  outside  the  orbit  of  the  earth.  This 
variation  in  parallactic  separation  yielded  heliocentric  distances  and  proper  inclinations 
for  the  observed  band  material  of  2.3  AU  and  8?7  (Gautier  et  al.  1984),  2.5  AU  and  8!l 
(Hauser  et  a/.1985),  and  2.44  AU  and  8°4  (Dermott  et  a/.1989).  None  give  error  estimates, 
but  the  values  seem  to  be  reasonably  consistent. 

Color  temperature  calculations  by  Low  et  al.  (1984)  produced  values  between  165 
K  and  200  K.  A  rapidly  rotating  graybody  of  this  temperature  would  be  located  between 
3.2  AU  and  2.2  AU,  well  within  the  main  asteroid  belt,  and  is  consistent  with  the  locations 
determined  by  parallax.  Consequently,  it  was  suggested  that  the  dust  bands  arose  from 
small  particles  generated  by  collisions  among  asteroids  (Low  et  al.  1984).  This  idea  was 
reinforced  by  an  apparent  association  between  the  latitudes  of  the  bands  and  the  proper 
inclinations  of  some  major  Hirayama  asteroid  families  (Dermott  et  a/.  1984),  which  are 
believed  to  have  been  produced  by  the  catastrophic  disruptions  of  large  asteroids  (Table 
II).  At  the  same  time,  Sykes  et  al. (1984)  calculated  that  the  random  catastrophic  disruption 
of  a  small  asteroid  (  10  km  in  diameter)  could  possibly  generate  sufficient  debris  to  be 

observed  by  IRAS  as  a  band  pair. 

Models  for  the  asteroidal  origin  of  the  dust  bands  has  since  followed  two  gen¬ 
eral  paths  which  will  be  examined  in  more  detail  later.  The  first  path  (Dermott  et 
a/.1984, 1985, 1986, 1988, 1989)  assumes  the  asteroid  belt  to  be  in  collisional  equilibrium, 
and  that  the  size  distribution  of  the  particle  population  at  all  locations  within  the  asteroid 


belt  is  characterized  by  a  power  law  with  a  single  index  over  all  sizes  from  tens  of  kilo¬ 
meters  J  t  tens  of  microns.  The  surface  area  of  dust  at  any  given  location  thus  increases 
monotonically  with  the  local  volume  of  asteroids.  Since  astere’d  families  represent  con¬ 
centrations  of  asteroids  in  a  —  sin  i  space,  there  should  be  corresponding  peaks  in  the  dust 
population.  Hence,  dust  bands  should  be  associated  with  asteroid  families. 

The  second  path  (Sykes  and  Greenberg  1986;  Sykes  1986,1988,1989)  is  a  non¬ 
equilibrium  theory  of  dust  band  origin  wh'ch  asserts  that  occasional  random  catastrophic 
disruptions  of  asteroids  results  in  debris  whose  subsequent  comminution  products  give  rise 
to  the  dust  bands.  In  this  case  dust  bands  are  not  associated  with  known  asteroid  families 
a  priori ,  and  may  be  found  at  other  locations.  This  model  predicts  that  dust  bands  should 
fade  with  time,  and  that  the  population  of  dust  bands  is  replenished  by  new  collisions. 
Though  focussing  on  the  more  frequent  disruptions  of  asteroids  tens  of  kilometers  in  di¬ 
ameter,  this  model  predicts  that  bands  arising  from  disruptions  large  enough  to  create  the 
largest  Hirayama  families  should  be  detectable  by  IRAS  for  >  1  Gyr. 

An  alternative  theory  for  dust  band  origin  is  the  disintegration  of  a  large  comet, 
which  will  be  considered  separately  in  this  chapter  (Sec.  III.C).  This  comes  about  from 
the  fact  that  the  potentially  most  prominent  cometary  suppliers  of  the  zodiacal  cloud  have 
similar  inclinations  to  dust  bands  (Dermott  et  a/.1984). 

In  the  next  section  a  mathematical  model  of  a  dust-band  torus  will  be  presented  to 
provide  some  insight  into  what  we  are  studying  in  terms  of  a  spatial  distribution  of  particles. 
The  effects  of  secular  gravitational  perturbations  and  dispersions  in  orbital  elements  on  this 
spatial  distribution  are  then  considered,  which  provide  additional  means  of  determining 
dust  band  locations  as  well  as  determining  the  effects  of  different  physical  processes  on  the 
small-particle  population  comprising  a  dust  band.  The  different  origin  scenarios  described 
above  are  presented,  along  with  their  observational  consequences.  Within  the  context  of 
the  (non-equilibrium)  random  collision  hypothesis,  we  examine  how  dust  bands  form,  how 
their  surface  areas  decrease  with  time,  and  how  many  dust  bands  are  likely  to  be  seen  as 
old  bands  are  replaced  by  newer  ones  from  more  recent  catastrophic  disruptions. 

II.  THE  DUST  BAND  TORUS 
A.  An  Analytical  Model 

The  spatial  distribution  of  the  dust  we  see  as  a  dust  band  pair  can  be  idealized  by 
considering  an  ensemble  of  dust  particles  whose  orbhs  have  identical  semi-major  axes  (a), 
proper  eccentricities  (e),  and  proper  inclinations  (t),  but  perihelia  and  nodes  distributed 
over  all  longitudes.  For  purposes  of  exposition  here,  we  neglect,  for  the  moment,  the  effects 
of  secular  gravitational  perturbations.  The  material  fills  a  torus  whose  particle  number 
density  is  (Sykes  1989) 
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within  the  limits 


a(l  —  e)  <  r  <  a(l  +  e)  (4) 

-i  <  /3  <  i  (5) 

Here  Cr  is  a  constant  of  normalization,  r  is  heliocentric  distance,  and  (3  is  the  heliocentric 
latitude  with  respect  to  the  plane  of  symmetry  (which  is  close  to  the  ecliptic  plane). 

These  equations  describe  a  spatial  distribution  that  is  a  torus  with  a  squarish 
cross-section  (Figure  3a),  centered  on  the  sun,  whose  radial  extent  is  bounded  by  the 
perihelion  and  aphelion  distances,  and  whose  latitudinal  extent  is  bounded  by  the  proper 
inclination  of  particle  orbits.  Maxima  in  volume  density  occur  at  the  “corners”  of  the  torus. 
The  latitudes  of  the  maxima,  seen  from  the  sun,  correspond  to  the  proper  inclination  of 
the  particle  orbits.  That  is  why,  from  the  earth,  these  concentrations  appear  as  pairs  of 
bands  straddling  the  ecliptic.  The  geometry  is  such  that  maxima  in  volume  density  at 
perihelion  and  aphelion  give  rise  to  both  perihelion  and  aphelion  band  pairs  (see  Permott 
et  al.  1985),  which  overlap  each  other  when  seen  from  the  Earth.  However,  IRAS  detects 
more  flux  from  perihelion  bands  for  two  reasons:  larger  particle  number  densities  and 
higher  temperatures. 

B.  The  Effects  of  Dispersions  in  Orbital  Elements 

Particles  that  make  up  a  dust  band  torus  do  not  have  identical  a,  e,  and  i;  there 
must  be  some  dispersion  in  orbital  elements.  This  significantly  affects  the  spatial  density 
of  dust  band  particles  (Figure  3b),  and  hence  the  locations  at  which  their  flux  is  seen. 
Two  mechanisms  by  which  the  orbital  elements  of  dust-band  particles  are  distributed  will 
be  examined.  The  first  mechanism  is  the  collisional  production  of  dust,  considered  in  the 
next  section.  The  second,  considered  in  the  section  following,  is  Poynting-Robertson  drag. 

Collisional  Dispersion.  Studies  of  collision  ejecta  indicate  that  smaller  particles  tend  to 
have  larger  ejection  velocities  (Melosh  1989),  which  yields  a  greater  dispersion  in  orbital 
elements  for  small  particles  than  large  particles.  Moreover,  in  a  plausible  particle  size 
distribution  (steeper  than  1/diameter),  small  particles  experience  more  collisions  with 
other  objects  of  comparable  size  than  do  large  particles.  The  orbital  elements  of  small 
particles  will  spread,  consequently,  more  rapidly  than  large  particles  The  most  obvious 
effects  are  associated  with  a  distribution  in  orbital  inclinations,  which  primarily  determines 
the  angular  separation  of  the  bands  we  observe.  For  a  Gaussian  distribution  with  fixed 
mean  values,  increasing  the  dispersion  in  proper  inclinations  results  in  (a)  increased  band 
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widths  and  (b)  shifting  peak  emission  of  the  latitudinal  profile  to  lower  latitudes  (Sykes 
1989;  Figure  4a).  Thus,  in  the  case  of  an  asteroid  family,  a  pair  of  bands  consisting  of 
associated  dust  having  submillimeter  diameters  would  be  expected  to  be  observed  at  lower 
latitudes  than  if  the  dust  band  particles  had  a  distribition  of  orbital  inclinations  identical 
to  kilometer  and  larger  family  members. 

In  general,  the  outer  edge  of  an  individual  band  (defined  as  the  latitude  where 
peak  flux  drops  by  half)  increases  slightly  (A  +  )  with  respect  to  the  mean  inclination  of 
particle  orbits  with  increasing  dispersion  in  inclination,  6i.  Similarly,  the  latitude  of  the 
peak  flux  and  inner  band  edge  (defined  as  the  latitude  where  the  flux  equals  the  average  of 
the  values  at  the  midplane  and  peak)  are  displaced  away  from  the  mean  inclination  towards 
the  midplane  by  Ao  and  A_,  respectively.  The  empirical  relations  describing  these  shifts 
are  (Sykes  1989): 
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Another  effect  that  can  be  seen  in  Figure  4  is  the  decreasing  contrast  of  the  inner  edges 
of  the  bands  with  increasing  dispersion  in  inclination.  The  ratio  of  peak  flux  ( Fp )  to  the 
flux  at  the  midplane  (Fc)  is  approximated  by 


Fp 

Fc 


(9) 


a  factor  of  2  lower  than  the  upper  limit  of  1.414/\/<5t'/i  determined  by  Dermott  et  al. (1985). 
The  relations  (6)  through  (9)  are  good  to  within  a  few  percent  for  i  <  20°  and  6i/i  <  0.25. 

Dispersions  in  semi-major  axis  and  eccentricity  act  to  increase  the  heliocentric 
distance  of  the  peak  flux  from  the  perihelion  bands  while  decreasing  it  for  the  aphelion 
bands  (Figures  4b  and  c).  Perihelion  and  aphelion  bands  associated  with  asteroid  families 
overlap  along  the  line  of  sight  when  observed  from  the  earth,  since  the  latitudinal  dirp'ace- 
ment  of  their  peaks  is  small.  Increasing  dispersion  in  semi-major  axis  and  eccentricity  thus 
decreases  an  already  small  peak-to-peak  separation. 

Dispersion  Due  to  Poynting- Robertson  Drag.  The  thermal  flux  from  the  dust  bands  is 
due  principally  to  those  particles  in  the  10  to  100  micron  size  range  because  of  their  greater 
surface  area  (Sykes  and  Greenberg  1986;  Dermott  et  al.  1986),  and  it  is  possible  that  the 
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orbits  of  these  particles  decay  significantly  due  to  Poynting-Robertson  drag  before  they  are 
either  dynamically  scattered  (Dermott  et  al.  1986)  or  are  comminuted  to  sizes  which  are 
ejected  from  the  immediate  vicinity  by  radiation  pressure  (Sykes  1989).  If  orbital  decay 
is  significiant,  then  the  particles  in  the  dust  bands  will  have  a  wide  range  of  semi-major 
axes,  and  may  show  significiant  dispersions  in  proper  eccentricity  and  inclination,  forced 
inclination,  and  forced  ascending  node  (these  “forced”  elements  are  functions  of  semi-major 
axis  and  will  be  discussed  more  fully  in  Sec.  II. C  below).  These  dispersions  arise  from  two 
causes:  (1)  Passage  through  various  resonances,  particularly  the  1:3  gap  at  2.5  AU,  results 
in  a  dispersion  of  the  proper  eccentricities  and  inclinations  (Dermott  et  a/.1989).  (2)  If 
the  particles  have  a  range  of  semi-major  axes,  then  they  must  also  have  a  range  of  forced 
orbital  elements  (cf.  Figure  5).  This  may  be  particularly  important  for  those  particles 
close  to  the  inner  edge  of  the  asteroid  belt  (a  <  2.5  AU),  and  it  may  be  that  the  dispersion 
of  the  orbital  elements  in  this  region  of  the  belt  defines  the  inner  edges  of  the  dust  bands 
if  the  particle  orbits  have  decayed  (Dermott  and  Nicholson  1989). 

The  effect  of  dispersion  due  to  Poynting-Robertson  drag  is  not  the  same  as  col- 
lisional  dispersion  (where  the  mean  orbital  elements  are  unchanged).  In  the  former  case, 
the  shift  in  the  peak  dust  band  emission  to  higher  or  lower  ecliptic  latitudes  is  a  function 
of  the  longitude  of  observation  (Sykes  1989).  Also,  a  simple  decrease  in  the  heliocentric 
distance  of  dust  percenters  results  in  an  increase  in  the  apparent  angular  separation  of  a 
band  pair  as  a  consequence  of  parallax. 

C.  The  Effects  of  Secular  Gravitational  Perturbations 

A  careful  examination  of  the  dust-band  images  in  ecliptic  coordinates  (Figure 
2)  shows  that  they  are  not  symmetric  about  the  plane  of  the  ecliptic.  The  orbits  of  the 
dust-band  particles  are  perturbed  by  Jupiter  and  other  planets  whose  orbits  are  slightly 
inclined  to  the  ecliptic.  Ensembles  of  particles  with  the  same  semi-major  axes  will  precess 
about  a  common  plane  where  the  torques  produced  by  these  gravitational  perturbations 
vanish.  This  defines  the  plane  of  symmetry  of  the  dust-band  torus  which  has  an  inclination 
relative  to  the  ecliptic)  i /  and  ascending  node  fl /  which  vary  with  semi-major  axis  (Figures 
5a  and  b).  These  gravitational  perturbations  also  act  to  distort  particle  orbits,  introducing 
a  “forced”  component  to  their  eccentricities.  Some  of  the  effects  of  this  on  the  dust-band 
torus  are  shown  in  Dermott  et  o/.(l985).  The  torus  shifts  away  from  its  sun-centered 
position  by  an  amount  ae/  in  a  direction  opposite  to  that  of  u)/  (Figure  5c),  which  now 
defines  the  longitude  of  pericenter  of  the  torus.  Figures  5b  and  d  show  the  variations  in 
these  parameters  as  a  function  of  semi-major  axis.  We  see  that  as  the  semi-major  axis 
moves  closer  to  Jupiter  (the  dominant  perturber)  beyond  2.6  AU,  these  forced  element 
components  approach  values  for  Jupiter’s  orbit.  The  shifting  of  the  torus  center  gives  rise 
to  a  longitudinal  temperature  variation  in  the  bands  which  are  diagnostic  of  their  distances 
(Dermott  et  al.  1985). 

The  existence  of  these  various  effects  allows  for  the  potential  extraction  of  all  the 
orbital  elements  of  the  dust  band  particles.  Semi-major  axis  can  be  determined  by  the 
direct  measurements  of  f!/  and  Uf.  The  use  of  fl/  in  this  way  is  described  for  the  (3  bands 
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in  Section  III. A.  The  apsidal  longitude  Cjj  is  that  at  which  dust-band  temperatures  are 
highest.  Given  the  semi-major  axis  of  the  torus  particle  orbits,  their  proper  eccentricity 
and  forced  eccentricity  can  be  decoupled  through  the  longitudinal  variation  of  dust  band 
temperature  (Dermott  and  Nicholson  1989).  The  proper  inclination  and  forced  inclination 
of  the  dust  bands  can  be  decoupled  from  measurements  of  band  center  latitudes  and 
band  separation  as  a  function  of  longitude  (Dermott  et  a/.1989;  Sykes  1989).  Secular 
gravitational  perturbations  also  result  in  the  distortion  of  the  dust  band  torus,  resulting 
in  the  north  and  south  bands  no  longer  being  exactly  plane  parallel.  Determining  the 
angle  between  the  planes  containing  the  individual  bands  making  up  a  band  pair  provides 
another  means  of  calculating  the  forced  eccentricity  (Sykes  1989). 

III.  THE  ORIGIN  OF  THE  DUST  BANDS 


A.  The  Collisional  Equilibrium  (Asteroid  Family)  Hypothesis 

This  hypothesis  assumes  that  the  dust  population  at  a  given  location  is  related 
to  the  population  of  observatble  asteroids  at  the  same  location.  This  has  been  examined 
by  assuming  that,  at  each  point  in  the  asteroid  belt,  all  sizes  of  particles  are  character¬ 
ized  by  a  single  equilibrium  size  distribution  arising  from  the  general  comminution  of  the 
asteroid  belt  through  mutual  collisions.  This  model  predicts  that  prominent  dust  bands 
are  associated  with  known  concentrations  of  asteroids  such  as  the  asteroid  families.  By 
demonstrating  a  relationship  between  dust  bands  and  asteroid  families,  the  equilibrium  hy¬ 
pothesis  provides  one  vehicle  by  which  dust  production  can  be  understood  in  the  asteroid 
belt  as  a  whole,  as  well  as  its  relationship  to  the  observed  zodiacal  dust  complex. 

The  relationship  between  the  Hirayama  asteroid  families  and  the  prominent  dust 
bands  has  been  a  central  question  in  the  study  of  the  latter  since  it  was  first  posed  by 
Dermott  et  ah(l984).  In  the  following,  the  consistency  between  the  calculated  surface  areas 
and  volumes  of  the  dust  bands  and  the  major  Hirayama  asteroid  families  is  examined.  Two 
basic  approaches  are  then  taken  to  determine  whether  the  dust  bands  derive  from  asteroid 
families.  The  first  method  is  the  direct  (or  indirect)  measurement  of  some  of  the  orbital 
elements  of  the  dust  bands.  These  can  then  be  compared  with  the  corresponding  elements 
of  the  asteroid  families.  The  second  method  is  the  generation  of  predictive  models  which 
are  then  compared  with  the  IRAS  dust  band  observations. 

Estimating  Dust  Band  Volume.  The  surface  optical  depth  of  the  prominent  dust  bands 
was  estimated  to  be  ~  10~8  (Low  et  ahl984),  corresponding  to  a  total  surface  area  of 
~  2  x  1019  cm2  (Dermott  et  ahl984).  It  is  assumed  that  the  size-frequency  distribution  of 
the  dust  band  particles  is  described  by  a  single  power  law  of  the  form 

d7V(m)  =  Km~qdm  (10) 

where  diV  is  the  number  of  particles  having  masses  between  m  and  m  +  dm,  K  is  a  constant, 
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and  q  is  the  mass  index.  Following  Dermott  and  Nicholson  (1989),  this  converts  to  the 
following  cumulative  power  law  in  radius, 


7V(r) 


__L _ fro 

3(q  -  1)  V  r 


3(9-1) 


(11) 


where  JV(r)  is  the  number  of  asteroids  with  radii  >  r  and  ro  is  a  constant,  similar  in 
magnituda  to  the  largest  particle  (in  this  case  asteroid)  radius.  The  total  area,  A,  of  the 
particles  is  then  given  (for  q  >  5/3)  by 


A  = 


7rr 
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(12) 


where  rmin  is  the  lower  cut-off  in  the  size  distribution,  while  the  total  volume  of  the 
particles,  V,  can  be  obtained  from  the  radius,  Re,  of  a  sphere  with  the  same  volume. 
V  =  (4/3)  7T Rg  and,  to  a  good  approximation  (for  q  <  2), 


R,  =  r0[3(2  -  «)]-,/3. 


(13) 


Values  of  Re  that  are  needed  to  account  for  the  observed  areas  are  shown  in  Figure  6  for  a 
range  of  values  of  rmtn  and  q.  The  theoretical  equilibrium  solution  for  asteroids  (Dohnanyi, 
1978)  has  a  population  index  q  =  1.837  (indicated  by  the  vertical  line  in  Figure  6).  If  q 
is  in  the  vicinity  of  this  value,  and  the  size-distribution  extends  down  to  ~  10^m,  then 
we  can  conclude  from  Figure  6  that  spheres  containing  the  inferred  volume  of  dust  band 
particles  (up  to  asteroids  tens  of  kilometers  in  radius)  must  have  equivalent  radii  of  ~  100 
km.  The  equivalent  radii  shown  in  Figure  6  are  similar  to  the  radii  of  Hirayama  family 
progenitors  (Gradie  et  al.  1979;  Dermott  et  al.  1984).  Assuming  a  single  size-frequency 
index,  however,  does  result  in  large  variations  in  surface  area  produced  by  a  given  mass 
of  material  if  the  index  or  minimum  particle  size  is  varied.  In  a  non-equilibrium  case 
(Section  III.B.;  Sykes  and  Greenberg  1986)  the  initial  population  index  may  be  close  to 
q  =  2.  Subsequent  evolution  calculated  by  Sykes  and  Greenberg  (1986)  shows  that  the 
population  quickly  changes  from  the  initial  power-law.  The  result  is  that  much  smaller 
bodies  (with  equivalent  radii  of  ~  10  km)  could  be  dust  band  parents,  and  numerical 
experiments  indicate  that  the  theory  is  not  very  sensitive  to  the  initial  value  of  q. 


Measuring  Dust  Band  Orbital  Elements.  The  simplest  means  of  determining  the  semi¬ 
major  axes  of  dust  band  particles  is  to  determine  the  orientation  and/or  the  inclination  of 
their  plane  of  symmetry  relative  to  the  ecliptic.  Utilizing  different  methods,  contradictory 
results  have  been  obtained  for  the  central  dust  bands. 

Using  the  0!5  IRAS  Zodiacal  History  File  (ZOHF)  (see  IRAS  Explanatory  Supp. 
1988),  Dermott  et  al. (1988)  separated  dust-band  profiles  from  the  broad  zodiacal  back¬ 
ground  utilizing  a  Fast  Fourier  Transform  and  high-pass  filtering  with  a  Parzen  window. 
Plotting  the  latitude  of  the  central  dust  band  (Figure  7),  in  which  the  a  and  (5  bands  are 
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not  resolved,  it  was  estimated  that  the  forced  inclination  of  the  composite  central  band  is 
1.2  and  that  the  longitude  of  the  associated  ascending  node  is  52!  The  forced  inclination 
clearly  differs  from  that  of  the  zodiacal  cloud  (l!5)  and  is  what  one  would  expect  for  the 
Themis  (at  a  =  3.1  AU)  and  Koronis  (at  a  =  2.9  AU)  asteroid  families  (see  Figure  5). 
However,  the  node  is  the  same  as  that  of  the  background  cloud  and  disagrees  with  the 
expected  longitude  of  ~  97°  for  particles  at  a  distance  of  3  AU  (Figure  5). 

A  parallactic  measurement  (Figure  8)  of  the  resolved  (3  band  pair  was  obtained 
by  Sykes  (1989)  using  the  IRAS  Skyflux  Maps  (IRAS  Explanatory  Supp.  1988).  The 
proper  inclination  of  the  particle  orbits  measured  was  2!l  ,  consistent  with  the  Koronis 
asteroid  family.  Sykes  found  the  ascending  node  for  this  band  pair  to  be  90°  ±  6°,  rea¬ 
sonably  consistent  with  the  value  of  96°  expected  for  particles  with  the  semi-major  axes 
of  the  Koronis  family.  The  discrepancy  with  the  earlier  measurement  is  thought  by  Sykes 
(1989)  to  arise  from  the  contamination  of  Fourier-filtered  dust-band  profiles  of  Dermott  et 
al.  (1988)  by  a  possible  cusp-like  high-spatial  frequency  component  of  the  broad  zodiacal 
emission. 

Comparison  with  Family  Dust- Band  Models.  Dermott  and  Nicholson  (1989)  have  con¬ 
structed  a  three-dimensional  numerical  model  that  permits  the  calculation  of  the  distri¬ 
bution  of  night-sky  brightness  that  would  be  produced  by  any  particular  distribution  of 
dust-particle  orbits.  This  model  includes  the  effects  of  planetary  perturbations  on  the  dust- 
particle  orbits,  reproduces  the  exact  viewing  geometry  of  the  IRAS  telescope,  and  allows 
for  the  eccentricity  of  the  earth’s  orbit.  The  result  is  a  model  for  the  variation  with  ecliptic 
latitude  of  the  brightness  observed  in  a  given  waveband  as  the  line  of  sight  of  the  telescope 
sweeps  through  the  model  dust  bands  at  a  constant  elongation  angle.  The  models  assume 
the  same  dispersion  in  orbital  elements  as  that  displayed  by  the  known  family  members. 
In  Figure  9  this  model  is  used  to  generate  a  predicted  profile  for  the  Eos  and  Themis 
asteroid  families  and  is  compared  with  a  high-pass  filtered  scan  from  the  Zodiacal  History 
File.  The  model  distributions  for  the  asteroid  families  are  too  sharply  peaked  to  account 
for  the  observed  profiles.  Dermott  and  Nicholson  then  examine  the  effects  of  dispersion 
due  to  evolution  under  Poynting- Robertson  drag,  and  this  yields  some  improvement  in  the 
comparison  of  model  and  filtered  profiles.  However,  while  the  agreement  between  observed 
dust  bands  and  specific  models  of  orbitally  evolved  particles  originating  in  the  prominent 
Hirayama  families  is  suggest'  *  .  definitive  results  from  this  approach  are  not  yet  available. 

Sykes  (1989)  ut;'i  •  the  dust  band  model  of  Sec.  II  and  maps  the  inner  and 
outer  edges  of  the  lalitui  "  :  profile  of  the  dust  bands  as  well  as  the  location  of  peak 
emission  onto  the  two-dimens:'.--  -.1  IRAS  Skyflux  Maps.  This  allows  the  a  and  /?  bands 
to  be  individually  studied  (E  i,  re  10a),  whereas  they  are  indistinguishable  at  the  lower 
resolution  of  the  Zodiacal  History  File.  Like  Dermott  and  Nicholson  (1989),  Sykes  initially 
assumes  that  the  mean  orbital  elements  of  the  dust  band  particles  are  the  same  as  the 
corresponding  family  members.  The  only  parameter  varied  is  the  dispersion  in  proper 
inclination  which  is  assumed  to  be  gaussian.  Separations  between  perihelion  and  aphelion 
bands  are  found  to  be  small,  and  increasing  dispersions  in  semi-major  axis  and  eccentricity 
only  makes  that  separation  smaller.  Mapping  routines  fully  reproduce  the  IRAS  pointing 
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geometry. 

Comparison  of  Koronis  and  Themis  models  with  the  (3  and  a  bands,  respectively, 
yields  excellent  fits  to  the  data  as  observed  in  the  IRAS  Skyftux  Maps  at  all  longitudes 
(Figure  10b, c).  This  assumes  increased  dispersions  in  proper  inclinations  of  less  than  a  fac¬ 
tor  of  two  for  the  model  Koronis  bands  and  a  little  greater  than  two  for  the  model  Themis 
bands.  Modifying  the  Koronis  and  Themis  family  models  by  decreasing  the  semi-major 
axis  and  other  associated  orbital  elements  (including  “forced”  components)  consistent  with 
orbital  decay  by  Poynting-Robertson  drag,  results  in  significant  divergence  from  the  ob¬ 
served  locations  of  the  a  and  (3  bands.  At  a  =  2.4  AU,  the  discrepancy  can  be  as  large 
as  l!5  in  latitude.  This  compares  with  a  width  of  ~  1/4°  for  the  (3  bands  in  the  IRAS 
Skyflux  plates. 

The  7  bands  are  very  closely  tracked  by  the  model  Eos  bands  in  the  filtered 
Zodiacal  History  File  images  (cf.  Figure  1),  indicating  an  orientation  and  inclination  of 
the  7  torus  that  may  be  similar  to  that  of  the  Eos  family.  However,  the  Eos  family  band 
model  yields  a  latitude  of  peak  emission  which  is  systematically  greater  than  that  observed 
for  the  7  bands.  Also,  the  pericenter  distances  of  the  bands  determined  by  parallax  (Section 
I.B)  indicate  a  value  ~  0.4  AU  smaller  than  that  of  Eos  family  members.  The  former  may 
be  resolved  by  increasing  the  dispersion  in  proper  inclinations  of  the  model  Eos  bands 
by  a  factor  of  ~  2.  This  also  results  in  significant  broadenning  of  the  bands  (Figure 
3) ,  which  seems  to  agree  with  the  several  degree  width  of  the  individual  7  bands  (Sykes 
1989).  Dispersion  due  to  Poynting-Robertson  drag  (section  II. B)  may  also  account  for 
band  broadenning  and  and  a  decrease  in  apparent  latitude  (Dermott  et  al.  1989),  though 
such  effects  are  viewing-geometry  dependent  and  may  not  occur  at  all  longitudes  observed 
(Sykes  1989).  On  the  other  hand,  the  apparent  discrepancy  in  pericenter  distances  may 
be  explained  by  the  Poynting-Robertson  drag.  An  alternative  explanation  in  both  cases  is 
that  the  7  bands  and  Eos  asteroid  family  are  not  causally  related. 

In  addition  to  the  more  prominent  a,  (3,  and  7  bands,  fainter  dust  bands  reported 
by  Sykes  (1988)  (Table  I)  were  compared  to  model  bands  associated  with  several  other 
known  asteroid  families  (Table  II).  This  resulted  in  a  possible  correlation  between  the  Io 
family  and  the  J/K  bands  (Sykes  1989).  However,  no  bands  were  detected  in  association 
with  the  Flora,  Nysa,  or  Maria  families.  Flora  and  possible  Maria  bands  had  been  predicted 
previously  on  the  basis  of  the  collisional  equilibrium  hypothesis  (Dermott  et  al.  1985).  Their 
absence  coupled  with  the  apparent  existence  of  bands  not  associated  with  any  prominent 
asteroid  families  present  a  challenge  to  the  collisional  equilibrium  hypothesis  of  dust  band 
origin. 

B.  The  Non-Equilibrium  (Random  Collision)  Hypothesis 

The  asteroid  belt  contains  tens  of  thousands  of  objects  in  mutual-crossing  orbits. 
On  average,  the  population  of  particles  may  be  described  by  some  equilibrium  size-frquency 
distribution  as  discussed  in  the  previous  section.  However,  when  two  asteroids  collide,  a 
small  area  of  orbital  element  space  will  be  filled  with  their  debris.  As  this  debris  experiences 
further  collisions  with  background  interplanetary  dust  particles  (IDP’s),  the  local  dust 


population  will  be  enhanced.  If  the  collision  is  large  enough,  the  resultant  dust  population 
deriving  from  the  comminution  of  its  debris,  may  have  enough  surface  area  to  be  detected 
by  an  instrument  such  as  IRAS.  Eventually,  removal  of  dust  through  radiation  forces  and 
the  continual  erosion  of  the  debris  population  results  in  the  decline  of  the  surface  area  of 
dust.  Then  at  another  location,  another  collision  takes  place  and  the  process  is  repeated. 

The  nonequilibrium,  or  random  collision,  hypothesis  of  dust  band  origin  views 
dust  bands  as  the  product  of  a  stochastic  process  in  which  they  are  created  and  destroyed 
over  geologic  time,  while  maintaining  a  certain  average  population  whose  spatial  distribu¬ 
tion  may  be  different  at  different  times.  This  is  in  contrast  to  the  equilibrium  hypothesis, 
in  which  the  dust  band  population  and  spatial  distribution  is  effectively  steady-state. 

In  the  following  subsections,  some  of  the  consequences  of  the  non-equilibrium 
model  are  examined:  band  pairs  must  form  on  finite  timescales;  their  surface  area  changes 
with  time  as  mass  is  collisionally  redistributed  from  larger  sizes  to  smaller,  a  pseudo¬ 
equilibrium  population  of  bands  must  be  maintained  due  to  constant  “gain”  and  “loss” 
rates.  Finally,  predictions  grounded  in  the  non-equilibrium  hypothesis  must  be  tested 
against  the  IRAS  data. 

Dust  Band  Formation.  The  collisional  disruption  of  an  asteroid  results  in  fragments 
having  a  dispersion  in  semimajor  axes  that  is  likely  small  compared  with  the  semimajor 
axis  of  the  parent  body.  For  example,  the  major  Hirayama  families  (Gradie  et  al.  1979) 
have  a  dispersion  of  semi-major  axes  corresponding  to  relative  orbital  velocities  of  ~  100 
m/s,  compared  with  absolute  orbital  velocities  of  ~  17  km/s.  In  fact,  it  is  a  mystery  how 
even  a  value  as  great  as  100  m/s  could  have  been  achieved.  Such  a  small  dispersion  in 
semimajor  axes  results  in  these  debris  being  distributed  around  the  orbit  of  the  parent 
body  on  time  scales  of  ~  103  years. 

The  orbits  of  asteroids  (and  therefore  debris  in  the  asteroid  belt)  experience  sec¬ 
ular  precession  of  their  apsides  and  nodes  as  a  consequence  of  gravitational  perturbations 
by  Jupiter  and  the  other  planets.  The  dust-band  torus  is  formed  as  the  orbits  of  collisional 
debris  precess  at  different  rates  due  to  small  differences  in  their  orbital  elements,  primar¬ 
ily  semi-major  axis,  so  that  with  time  their  nodes  spread  around  the  ecliptic  relative  to 
each  other.  Two  mechanisms  operate  simultaneously  to  distribute  the  nodes  of  particle 
orbits.  The  first  is  differential  precession  due  to  the  differences  in  semi-major  axes  for 
particles  of  a  given  size.  The  second  is  a  mass  fractionation  of  the  nodes  arising  from  the 
size-dependent  variation  of  semi-major  axes  with  time  due  to  Poynting- Robertson  drag  on 
smaller  particles  (Sykes  and  Greenberg  1986). 

Increasing  the  dispersion  in  semimajor  axes  increases  the  rate  at  which  band  pairs 
are  formed  by  this  mechanism,  as  does  increasing  the  mean  semimajor  axis  of  the  particle 
orbits.  This  last  is  a  consequence  of  the  stronger  gravitational  effects  of  Jupiter  as  the 
particle  orbits  approach  Jupiter’s  orbit.  From  Figure  11,  minimum  formation  timescales 
range  between  105  and  ~  107  years,  with  dust  bands  associated  with  the  Themis,  Koronis, 
and  Eos  families  forming  in  less  than  106  years. 

Mass  fractionation  of  the  orbital  nodes  occurs  as  Poynting-Robertson  drag  de¬ 
creases  the  semi-major  axis  of  small  particles  with  time,  relative  to  large  particles.  Small 
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particles  consequently  precess  at  increasingly  lower  rates  with  respect  to  large  particles. 
Depending  on  the  size-frequency  distribution  of  particles  created  in  the  initial  catastrophic 
disruption,  this  mechanism  can  potentially  result  in  more  rapid  band  pair  formation  (~  105 
to  106  years)  than  differential  precession  due  to  the  initial  dispersion  in  semi-major  axes 
(Figure  12). 

Surface  Area  Evolution.  The  debris  that  comprises  a  dust-band  torus  will  not  maintain 
its  original  mass  distribution.  Particles  are  comminuted  by  collisions  with  background 
interplanetary  dust  particles,  both  cometary  and  asteroidal  in  origin.  Catastrophic  frag¬ 
mentation,  rather  than  gradual  erosion,  dominates  the  comminution  (Dohnanyi  1978), 
and  the  breakup  of  larger  particles  into  smaller  particles  increases  the  surface  area  of  the 
dust  bands,  tending  to  increase  their  brightness.  Mass  (and  surface  area)  from  the  torus 
is  eventually  lost  when  the  debris  is  comminuted  into  small  enough  fragments  that  they 
are  rapidly  swept  away  by  Poynting-Robertson  drag  or  radiation  pressure  (Figure  13). 
For  most  meteoritic  compositions  this  corresponds  to  particle  diameters  between  approxi¬ 
mately  0.08  and  4  microns  (Burns  et  al.  1979). 

In  the  competition  between  production  and  loss  rates  of  surface  area  in  a  dust 
band  torus,  gradual  loss  eventually  wins  out  with  time  (Figure  13c).  Thus,  the  eventual 
fate  of  any  band  pair  in  the  non-equilibrium  model  is  to  fade  away  (Sykes  and  Greenberg 
1986). 

How  Many  Dust  Bands  Should  We  See?  According  to  the  nonequilibrium  model,  the 
population  of  bands  above  a  limiting  surface  area  will  be  determined  by  two  quantities: 
the  mean  time  between  collisions  occur  which  generate  more  than  the  minimum  required 
surface  area  (the  gain  term),  and  the  lifetime  of  bands  before  they  fade  below  this  limit 
(the  loss  term) .  The  ratio  of  the  gain  to  loss  terms  determines  the  number  of  bands  one 
would  expect  to  observe  at  any  given  time.  Based  on  their  models,  Sykes  and  Greenberg 
(1986)  calculated  that  there  should  be  on  the  order  of  2  pairs  of  bands  with  surface  areas 
equal  to  or  exceeding  that  estimated  for  the  bands  reported  by  Low  et  al.  (1984)  and  ~30 
pairs  of  bands  an  order  of  magnitude  fainter.  This  agreed  nicely  with  the  number  of  bands 
at  the  bright  end  (Low  et  a/.1984),  but  seemed  to  conflict  with  the  smaller  number  of  faint 
bands  detected  by  Sykes  (1988).  The  problem  was  resolved  when  it  was  found  that  the 
dust  bands  are  not  sharp  features,  but  have  typical  widths  of  a  few  degrees  (the  exception 
being  the  a  and  f3  bands.  Assuming  Gaussian  profiles,  Sykes  (1988)  showed  that  the 
superposition  of  so  many  faint  and  bright  bands  would  produce  something  of  a  continuum 
above  which  only  2  to  5  bands  (or  clusters  of  bands)  could  be  detected  in  addition  to  the 
a,  (3,  and  7  bands.  Thus,  the  noneqilibrium  model  is  consistent  with  the  numbers  of  bands 
now  detected. 

The  dust  bands  predicted  by  Sykes  and  Greenberg  (1986)  were  disruption  prod¬ 
ucts  of  asteroids  tending  to  have  diameters  between  5  and  10  km.  A  significant  fraction 
of  these  bands  also  arose  from  the  disruption  of  30  to  90  km  diameter  bodies.  Thus,  the 
prominent  a,  /?,  and  7  bands  (with  the  possible  exception  of  the  /?  bands),  assuming  as¬ 
sociation  with  the  Themis,  Koronis,  and  Eos  families,  respectively,  were  not  predicted  by 
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Sykes  and  Greenberg  (1986).  Bodies  large  enough  to  form  the  Themis  and  Eos  asteroid 
families  disrupt  so  infrequently  (on  time  scales  of  the  age  of  the  solar  system)  that  they 
were  not  included  in  the  collision  statistics.  However,  the  Sykes  and  Greenberg  model  can 
explain  the  observation  of  these  bands  in  the  event  of  such  large  collisional  disruptions. 
This  is  due  to  their  long  observational  lifetimes.  At  the  other  end  of  the  scale,  asteroids 
much  smaller  than  5  km  in  diameter  would  not  produce  enough  debris  to  form  a  detectable 
band  pair. 

The  mean  time  between  catastrophic  collisions  generating  a  mass  of  debris  equiv¬ 
alent  to  a  5-10  km  asteroid  was  found  to  be  10s  and  106  years  -  on  the  order  of  both  the 
dust  band  formation  time  and  lifetime  (for  surface  areas  exceeding  1018  cm2,  ~  1/10  that 
of  the  prominent  bands).  Consequently,  it  was  thought  that  faint  partial  band  pairs  might 
be  observed  (whose  nodes  were  not  distributed  over  all  latitudes).  Such  a  partial  band  was 
reported  to  be  detected  (Sykes  1986,  1987),  but  later  evidence  suggested  that  the  structure 
had  moved  over  a  period  of  eight  months,  indicating  that  it  may  actually  be  a  debris  trail 
of  a  type  not  previously  observed  (Sykes  1988). 

C.  The  Comet  Hypothesis 

It  cannot  be  concluded  a  priori  that  because  the  zodiacal  dust  bands  are  located 
in  the  asteroid  belt,  that  they  have  an  asteroidal  origin.  Numerous  short-period  comets 
have  low  inclinations;  the  peak  in  the  distribution  of  the  observed  inclinations  is  actually 
close  to  10°.  Many  of  the  observed  short-period  comets  also  have  perihelion  distances 
between  2  and  4  AU  and,  because  of  selection  effects,  there  is  probably  a  large  number 
of  such  comets  with  perihelia  within  the  asteroid  belt  that  remain  undetected  (Burns  et 
al.  1984;  Dermott  et  al.  1989). 

IRAS  discovered  that  n.  y  short-period  comets  have  associated  dust  trails  de¬ 
tected  over  large  portions  of  their  orbits  (Sykes  et  al.  1986a).  In  Figure  1,  the  dust  trails 
associated  with  P/Tempel  2  and  P/Encke  are  clearly  seen.  This  represents  a  possibly 
significant  increase  in  the  contribution  of  cometary  dust  to  the  zodiacal  dust  complex, 
particularly  in  the  submillimeter  and  millimeter  size  ranges  (Sykes  et  al.  1986b).  Another 
source  of  zodiacal  dust  of  possible  cometary  origin  are  the  Type  II  dust  trails  (A  to  D 
in  Figure  1)  (Sykes  1988).  These  recently  detected  structures  may  arise  from  the  ejection 
of  large  particles  at  high  velocities  (~  100  m/s)  from  a  comet  nucleus  during  perihelion 
passage,  in  comparison  to  the  lower  velocities  (<  10  m/s)  associated  with  the  originally 
detected  Type  I  trails.  They  may  also  represent  the  breakup  of  cometary  nuclei.  In  either 
case,  the  dust  trail  particles  will  tend  to  evolve  into  a  dust  band  torus  as  described  in  Sec. 
III.B.,  though  with  a  much  greater  radial  width  than  for  asteroid  debris  due  to  the  larger 
typical  cometary  eccentricities. 

As  a  torus  forms,  the  cometary  material  is  spread  out  over  a  substantially  larger 
volume,  tending  to  decrease  the  surface  brightness  of  the  dust  observed  unless  it  is  increased 
through  comminution  of  the  larger  particles  or  the  emission  of  additional  material  from  the 
parent  comet.  In  the  case  of  the  latter,  emissions  from  short-period  comets  are  expected 
to  last  only  about  10,000  orbits  (the  nominal  lifetime  of  a  short-period  comet),  which  is 
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smaller  than  the  dust  band  formation  timescales.  Unless  the  comet  itself  is  disrupted,  the 
population  of  large  particles  whose  comminution  provides  the  observed  dust  population  is 
likely  to  be  minimal.  Therefore  dust  bands  associated  with  cometary  emission  are  likely 
to  be  very  short-lived.  There  remains  the  possibility,  however,  that  emissions  from  a  large 
number  of  small,  undetected  short-period  comets  could  superpose  to  yield  an  observable 
dust  band  torus. 

The  disruption  of  a  cometarj^  nucleus  might  allow  for  a  dust  band  to  be  formed 
as  bright  as  any  detected  by  IRAS.  Though  highly  uncertain,  the  size-frequency  distri¬ 
bution  of  the  large-particle  population  of  comets  is  thought  to  have  a  size  index  q  ~  2.1 
(Sekanina  1979).  This  is  similar  to  the  initial  size  distribution  assumed  for  a  catastroph¬ 
ically  disrupted  asteroid  by  Sykes  and  Greenberg  (1986).  From  Figure  6,  such  breakup 
would  require  an  initial  cometary  radius  of  only  ~  10  km,  similar  to  that  observed  for 
Halley  (Sagdeev  et  a/.1986;  Keller  et  a/.1986).  The  frequency  with  w'hich  such  cometary 
disruptions  occur  would  then  determine  the  probability  of  the  resultant  band  pair  being 
observed. 

One  means  of  distinguishing  between  asteroidal  and  cometary  models  of  the 
bands  is  probably  a  combination  of  the  forced  orbital  element  and  parallax  methods  (Der- 
mott  et  a/.1989).  The  forced  orbital  element  method  determines  the  semi-major  axis  of 
the  particles,  whereas  the  parallax  method  determines  an  “effective  distance”  which  is 
probably  closely  related  to  the  percenter  distance,  a(l  -  e).  Thus,  it  should  be  possible,  in 
principal,  to  determine  the  proper  eccentricities  of  the  particle  orbits,  thereby  indicating 
whether  the  particles  derive  from  comets. 

IV.  DISCUSSION 


Neither  the  equilibrium  nor  non-equilibrium  models  of  dust  band  origin,  as  cur¬ 
rently  formulated,  present  a  complete  picture  of  the  IRAS  dust  band  observations.  The 
equilibrium  model  fails  to  explain  the  existence  of  bands  at  non-family  locations  and  the 
absence  of  bands  associated  with  the  Flora,  Nysa,  and  Maria  families.  Both  theories  ex¬ 
plain  the  a,  /?,  and  7  bands  and  their  respective  relationships  to  the  Themis,  Koronis,  and 
Eos  families.  The  non-equilibrium  model  of  Sykes  and  Greenberg  (1986)  overestimates  the 
surface  brightness  of  bands  arising  from  the  disruption  products  of  small  asteroids  (un¬ 
less  it  eventually  turns  out  that  7  bands  are  not  associated  with  the  Eos  family).  It  also 
predicts  the  detection  of  partial  bands,  which  have  yet  to  be  unambiguously  observed. 

On  the  other  hand,  the  equilibrium  model  led  to  the  initial  correlation  between 
the  prominent  bands  and  the  major  Hirayama  asteroid  families  (Dermott  et  a/.1984),  and 
the  non-equilibrium  model  inspired  the  search  for  and  detection  of  zodiacal  structures 
including  additional  dust  bands  (Sykes  1988).  Both  models  predicted  that  a  significant 
fraction  of  the  zodiacal  dust  complex  derives  from  the  asteroid  belt. 

If  the  asteroid  belt  is  a  principal  source  of  zodiacal  dust,  then  that  dust  must 
be  transported  to  the  inner  solar  system  by  Poynting-Robertson  drag.  Initially,  it  seems 
contradictory  that  the  previous  statement  be  true  while  the  dust  bands  -  regions  where 
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dust  is  obviously  being  generated  -  evidence  little  if  any  orbital  decay  (Hauser  et  al.  1985; 
Dermott  and  Nicholson  1989;  Sykes  1989).  Sykes  (1989)  interprets  these  observations 
as  indicating  particle  sizes  whose  collisional  lifetimes  are  short  in  comparison  with  their 
Poynting- Robertson  decay  times,  and  that  the  removal  of  their  comminution  products 
from  the  ensemble  of  dust  band  particles  must  be  rapid.  This  does  not  mean  that  the 
particles  simply  disappear,  rather  that  they  are  no  longer  distinguishable  from  the  broad 
zodiacal  background  depicted  in  Figure  1.  One  possibility  is  that  a  significant  fraction 
of  the  smaller  particles  are  comminuted  into  sizes  sensitive  to  radiation  pressure,  which 
are  then  essentially  removed  from  the  ensemble  of  dust-band  particles  on  timescales  of 
an  orbital  period.  These  particles  eventually  decay  by  P-R  drag,  but  their  surface  area 
is  spread  over  a  large  volume  and  hence  contribute  little  to  the  local  dust-band  surface 
brightness. 

Dermott  and  Nicholson  (1989)  believe  that  P-R  drag  will  turn  out  to  explain  the 
broadness  of  the  bands  and  the  apparent  displacement  of  parallactic  distances  of  the  dust 
bands  inward  from  the  pericenters  of  their  associated  families. 

The  association  of  dust  bands  with  asteroid  families  may  help  to  shed  some 
light  on  the  nature  and  evolution  of  the  latter.  The  major  Hirayama  families  are  thought 
to  have  originated  with  the  breakup  of  large  parent  asteroids  (Kuiper  1950).  That  the 
dust  bands  have  orbital  elements  consistent  with  the  mean  elements  of  the  corresponding 
families  indicates  (from  the  perspective  of  the  non-equilibrium  model)  that  we  probably  are 
looking  at  comminution  products  of  debris  from  the  original  family-forming  collision,  rather 
than  from  the  random  disruption  of  a  smaller  family  member  (which  would  likely  have  a 
different  inclination).  Thus,  the  dust-band  surface  area  evolution  model  can  be  utilized 
in  an  attempt  to  constrain  the  ages  of  asteroid  families.  This  was  done  by  Sykes  (1986) 
utilizing  the  dust  band  surface  area  evolution  model  of  Sykes  and  Greenberg  (1986).  A 
total  surface  area  of  2  x  1019  was  assumed  for  bands  associated  with  the  Themis,  Koronis, 
and  Eos  families,  which  is  probably  something  of  an  overestimate,  particularly  for  the 
Themis  family  which  have  a  relatively  low  surface  brightness  in  the  IRAS  skyflux  maps. 
A  surface  area  of  10 18  cm2  was  assumed  as  an  upper  limit  to  the  dust  in  the  Nysa  family. 
The  calculated  times  elapsed  since  the  disruption  of  the  family  parent  bodies  are  listed 
in  Table  III,  and  should  be  considered  as  very  model  dependent.  If  the  calculated  values 
are  correct,  this  model  indicates  that  the  Koronis  family  is  much  younger  than  both  the 
Themis  and  Eos  family.  This  result  is  consistent  with  rotational  studies  of  family  members 
(Binzel  1987).  If  the  Nysa  family  was  the  product  of  a  catastrophic  disruption,  then  its 
lack  of  an  associated  dust  band  suggests  that  a  family-forming  collision  was  indeed  ancient. 

The  Sykes  and  Greenberg  model  says  that  dust  bands  deriving  from  the  dis¬ 
ruption  of  an  asteroid  on  the  order  of  a  few  hundred  kilometers  in  diameter  should  be 
detectable  by  a  detector  like  IRAS  for  at  least  1  Gyr.  Also,  according  to  this  model  the 
probability  of  any  such  event  occuring  in  the  last  several  Gyr  is  very  small.  The  fact  that 
two  such  events  are  observed  (Themis  and  Eos)  and  that  other  families  with  similar  sized 
parents  arose  from  even  earlier  collisions  suggests  two  things:  (1)  collisional  activity  was 
much  greater  in  the  asteroid  belt  more  than  a  billion  years  ago,  or  (2)  many  families  may 
not  have  derived  from  a  single  collisional  event,  but  that  they  may  be  associations  aris- 
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ing  from  segregation  or  lumping  (in  orbital  element  space)  as  a  consequence  of  dynamical 
resonances,  such  as  the  Phocaea  family  (Williams  1971). 

IRAS  first  detected  the  dust  bands  in  1983,  providing  a  new  database  and  phe¬ 
nomenology  against  which  models  of  asteroid  collisional  activity  and  small  particle  dy¬ 
namics  could  be  tested.  A  wealth  of  more  information  is  yet  to  come  as  more  spacebased 
infrared  telescopes  are  launched  (beginning  with  the  Cosmic  Background  Explorer  (COBE) 
in  the  fall  of  1989,  to  be  followed  by  the  Infrared  Space  Observatory  (ISO)  in  the  1990s 
and  the  Space  Infrared  Telescope  Facility  (SIRTF)  early  in  the  next  century).  The  instru¬ 
ments  being  launched  will  allow  not  only  for  more  detailed  studies  along  lines  reviewed 
here,  but  also  will  allow  spectroscopy  studies  to  be  undertaken,  opening  whole  new  areas 
of  investigation  to  link  dust  to  asteroids  and  comets  and  link  dust  in  the  asteroid  belt  to 
extraterrestrial  dust  collected  in  the  Earth’s  atmosphere  and  elsewhere. 
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Table  I.  Dust  Band  Pairs 


Name 

Geocentric 

Comments 

Ecliptic  Latitude 
(HCON  1  and  2) 

“ 

a 

±(0  to  2.5) 

Unresolved  with  0  (Low  et  al.  1984).  Resolved  ans  associated  with 
Themis  family  (Sykes  1986). 

0 

±(l  to  3.5) 

Unresolved  with  a  (Low  et  a/.1984).  Resolved  and  associated  with 
Koronis  family  (Sykes  1986). 

1 

±(8.5  to  11.5) 

Low  et  al.,  1984.  Associated  with  Eos  family  (Dermott  et  a/.1984). 

E 

4  to  6 

Possible  pair  with  F  (Sykes  1988) 

F 

-4  to  -6.5 

Possible  pair  with  E  (Sykes  1988) 

G 

6.5  to  8 

Possible  pair  with  H  (Sykes  1988). 

H 

-5.5  to  -8.5 

Possible  pair  with  G  (Sykes  1988). 

J 

12.5  to  15 

Probable  pair  with  K  (Sykes  1988).  Associated  with  Io  family  (Sykes 
1989). 

K 

-13  to  -16 

Probable  pair  with  J  (Sykes  1988).  Associated  with  Io  family  (Sykes 
1989). 

M 

15  to  17.5 

Probable  pair  with  N  (Sykes  1988). 

N 

-17  to  -20 

Probable  pair  with  M  (Sykes  1988). 

Table  II.  Asteroid  Family  Proper  and  Forced  Orbital  Elements 


Family 

a 

6a 

e 

6e 

ef 

i 

6i 

lf 

nf 

Themis 

3.137 

0.025 

0.153 

0.009 

0.038 

1.420 

0.223 

1.22 

97. S 

8.7 

Koronis 

2.875 

0.018 

0.049 

0.006 

0.037 

2.118 

0.081 

1.16 

96.1 

G.2 

Nysa 

2.446 

0.020 

0.162 

0.009 

0.036 

3.205 

0.154 

0.93 

86.9 

351.7 

Flora 

2.193 

0.008 

0.138 

0.006 

0.048 

5.024 

0.452 

0.61 

43.4 

338.2 

Eos 

3.015 

0.006 

0.071 

0.008 

0.037 

10.12 

0.710 

1.19 

97.1 

7.6 

Io 

2.650 

0.035 

0.143 

0.009 

~  0.06 

13.36 

0.240 

1.06 

93.9 

~  32 

Maria 

2.550 

0.013 

0.089 

0.009 

0.035 

15.21 

0.170 

1.03 

91.0 

354.2 
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Table  III.  Asteroid  Family  Formation  Estimates 


Family 

Minimum 
Parent  Body 
Diameter 
(km) 

Minimum 

Age 

(years) 

Themis 

300 

1.5  x  109 

Koronis 

90 

1.3  x  108 

Nysa 

200 

2  x  109 

Eos 

189 

5.8  x  108 
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Figure  Captions 


Figure  1.  Brightness  profiles  for  an  IRAS  scan  on  1983  June  24  at  elongation  91°  1  (Hauser 
et  al,  1984).  The  ecliptic  plane  is  crossed  at  longitude  l!3.  The  galactic  plane  is  crossed 
near  ecliptic  latitude  +60°  and  galactic  longitude  l  =  96!  8.  All  discrete  features  are  from 
real  sources.  The  prominent  source  near  the  south  ecliptic  pole  at  60  and  100  //m  is  the 
Large  Magellanic  Cloud.  At  12  /zra  and  25  the  thermal  emission  is  dominated  by  a 
broad  zodiacal  component  while  galactic  emission  dominates  at  the  longest  wavelengths. 
The  arrows  indicate  bumps  in  the  profile  corresponding  to  the  -7  bands. 

Figure  2.  The  25/xm  IRAS  scans  have  been  segregated  into  (a)  HCON  2  and  (b)  HCON  3. 
and  zero-sum  high-pass  filtered  along  the  in-scan  direction  using  a  filter  width  of  3! 5.  The 
data  has  been  binned  into  0!5  x  0!5  pixels  and  mapped  in  ecliptic  cylindrical  coordinates. 
Because  of  the  separation  of  only  ~  10  days,  HCON  1  (not  shown)  and  HCON  2  images 
are  nearly  identical.  The  dust  band  structures  are  described  in  Table  1.  A  —  D  are  Type 
II  dust  trails,  T2DT  and  EDT  denote  the  Tempel  2  and  Encke  dust  trails,  respectively, 
and  GC  indicates  galactic  cirrus.  The  galactic  plane  is  seen  crossing  the  ecliptic  near  90° 
and  270°  ecliptic  longitude  (Sykes  1988). 

Figure  3.  (a)  The  particle  number  density  distribution  in  a  radial  cross  section  of  a 
Koronis  dust  band  torus  whose  particles  have  orbital  elements  equal  to  the  mean  values 
of  the  Koronis  asteroid  family  (with  no  dispersion  in  elements).  The  Sun  is  to  the  left. 
The  jagged  edges  are  an  artifact  of  sampling,  (b)  Convolution  by  Gaussian  dispersions  in 
proper  inclination  (la)  and  semi-major  axis  (2a)  of  the  model  Koronis  torus  increases  the 
volume  of  the  torus,  particularly  in  the  radial  direction,  while  separations  between  peaks 
densities  at  perihelion  and  aphelion  decrease  both  radially  and  in  latitude  (Sykes  1989). 

Figure  4.  (a)  The  latitudinal  profile  of  the  dust  bands  are  shown  after  being  convolved 
with  a  gaussian  distribution  of  orbital  inclinations  with  standard  deviation  Si.  Increasing 
dispersion  in  inclination  results  in  decreasing  latitudinal  separation  of  the  dust  bands  and 
decreasing  contrast  between  the  band  peaks  and  centers.  Ao,  A  +  ,  and  A_  indicate  the 
shift  of  peak  and  band  edge  latitudes,  respectively,  and  are  defined  in  the  text.  The  mean 
elements  of  the  particles  remain  unchanged,  (b)  Radial  profiles  of  the  particle  number 
density  of  a  model  torus  associated  with  the  Koronis  asteroid  family  are  shown  for  different 
dispersions  in  semimajor  axis  and  (c)  eccentricity.  Densities  are  normalized  to  the  value 
at  a  =  2.875  AU.  The  la  values  correspond  to  the  standard  deviation  of  those  terms  over 
the  orbits  of  known  family  members  (Table  II,  calculated  from  Williams  1979).  Increasing 
dispersion  in  both  terms  results  in  a  decrease  in  the  radial  separation  of  the  perihelion  and 
aphelion  bands,  as  well  as  a  decrease  in  their  contrast.  The  mean  elements  of  the  particles 
are  unchanged  (Sykes  1989). 

Figure  5.  (a)  Variation  of  the  ecliptic  latitude  of  a  dust  band  with  ecliptic  longitude,  as 
measured  in  a  sun-centered  coordinate  system.  The  latitudinal  width  of  the  band  at  all 
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longitudes  is  2 i',  twice  the  proper  inclination  of  the  dust  particle  orbits.  The  amplitude  of 
the  apparent  sinusoidal  variation  is  determined  by  the  forced  inclination,  v j.  (b)  Variation 
of  the  forced  inclination  if  and  the  corresponding  longitude  of  the  ascending  node  fl /  with 
semi-major  axis.  The  elements  with  subscript  J  are  the  present  orbital  elements  of  Jupiter, 
(c)  Plane  view  of  the  distribution  of  elliptical  particle  orbits  in  a  dust  band.  To  order  e, 
the  figure  is  circularly  symmetric  about  the  point  C  which  is  displaced  from  the  sun  by  an 
amount  atf  in  a  direction  opposite  to  that  of  the  forced  pericenter,  cD/.  (d)  The  variation 
of  the  forced  eccentricity  e/  and  the  corresponding  longitude  of  the  forced  pericenter,  u )/, 
with  semi-major  axis  (Dermott  et  al.,  1985). 

Figure  6.  The  total  surface  area  of  particles  contained  within  a  sphere  of  equivalent 
radius,  Re,  is  shown  for  a  range  of  power  law  distribution  index,  q,  and  the  lower  cut-off 
in  the  particle  size  distribution,  rmin.  The  horizontal  line  shows  the  area  for  the  central 
dust  band,  inferred  from  Low  et  al. (1984)  (Dermott  and  Nicholson  1989). 

Figure  7  Variation  of  the  latitude  of  the  center  (measured  between  the  half-power  points) 
of  the  central  dust  band  as  with  day  of  the  year,  as  seen  at  25  yum  in  the  IRAS  Zodiacal 
History  File.  Data  obtained  when  IRAS  was  in  the  leading  or  ascending  leg  of  its  motion 
around  the  earth  are  represented  by  squares,  while  the  data  from  the  trailing  or  descending 
leg  are  represented  by  triangles.  The  vertical  lines  give  the  forced  inclination  of  the  central 
band  and  the  location  of  the  nodes  (Dermott  et  al.  1988). 

Figure  8  Geocentric  ecliptic  coordinates  of  the  north  and  south  (3  bands  as  measured  from 
unfiltered  IRAS  skyflux  plates  (25  ptm).  The  latitude  of  each  band-pair  component  was 
binned  and  averaged  every  0!5  in  longitude  (filled  circles).  North  and  South  bin  values 
were  then  averaged  (open  circles)  in  order  to  determine  the  geocentric  ecliptic  longitude 
of  the  ascending  node  as  seen  in  both  HCON  1  and  HCON  3,  roughly  eight  months  apart. 
The  earth’s  longitude,  A®,  is  given  for  the  times  of  plane-crossing  observations.  HCON  1 
observations  plotted  were  made  in  the  direction  opposite  to  the  earth’s  motion,  while  the 
HCON  3  observations  plotted  were  made  in  the  direction  of  earth’s  motion  (Sykes  1989). 

Figure  9  An  observed  dust  band  profile  (smooth  curve),  after  high-pass  filtering,  is  com¬ 
pared  with  the  prediction  of  a  model  based  on  the  distribution  of  orbits  in  the  Fos  and 
Themis  families  (Dermott  and  Nicholson  1988). 

Figure  10  (a)  An  IRAS  skyflux  map  (Plate  95,  25/zm)  has  been  boxcar  high-pass  filtered, 
using  a  filter  width  of  1°  along  the  scan  direction,  allowing  the  a  and  f3  bands  to  be  dis¬ 
tinguished.  The  edges  of  model  (b)  Koronis  and  (c)  Themis  dust  bands  are  then  projected 
onto  the  skyflux  map  for  comparison  with  the  observed  bands  (Sykes  1989). 

Figure  11.  The  longitudes  of  the  ascending  nodes  of  an  ensemble  of  collisional  debris 
orbits  are  increasingly  dispersed  due  to  differences  in  semi-major  axes  among  the  debris 
after  ejection.  Nodal  dispersion  as  a  function  of  time  is  shown  for  different  mean  ejection 
velocities  in  the  inner  (a)  and  outer  (b)  part  of  the  asteroid  belt.  When  the  nodes  have 
been  distributed  over  360°,  band  formation  is  complete.  Band  formation  is  more  rapid  in 
the  outer  belt  due  to  the  closer  proximity  of  Jupiter  (Sykes  and  Greenberg  1986). 
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Figure  12.  The  ascending  nodes  of  debris  orbits  from  a  single  collisional  event  are  shown 
as  a  function  of  time.  Poynting-Robertson  drag  coupled  with  Jovian  gravitational  pertur¬ 
bations  acts  to  decrease  the  rate  of  orbital  precession  by  small  particles  relative  to  large 
particles  by  decreasing  their  semi-major  axes  with  time.  This  results  in  band  pair  forma¬ 
tion  by  mass  fractionation  of  debris  orbits.  The  distribution  of  ascending  nodes  by  this 
mechanism  (thin  curves)  and  by  the  initial  dispersion  of  semi-major  axes  (thick  curves) 
are  shown  at  different  times  for  a  single  collision  in  the  inner  (a)  and  outer  (b)  part  of 
the  asteroid  belt.  Particle  mass  densities  of  p  =  3  g/cm3  are  assumed,  and  a  dispersion  in 
ejection  velocities  of  100  m/s.  For  reference,  the  inner  circle  represents  the  orbit  of  earth. 
As  in  Figure  11,  band  pair  formation  in  the  outer  belt  is  more  rapid  than  the  inner  belt 
due  to  the  closer  proximity  of  Jupiter  (Sykes  and  Greenberg  1986). 

Figure  13(a)  Debris  from  the  catastrophic  disruption  of  an  asteroid  has  an  initial  size- 
frequency  distribution  which  subsequently  (b)  undergoes  collisional  redistribution  of  its 
mass  from  larger  to  smaller  particles  until,  eventually,  particles  are  small  enough  that  they 
are  rapidly  removed  by  radiation  forces  from  the  ensemble  comprising  a  dust  band  torus, 
(c)  The  surface  area  of  dust  band  tori  gradually  diminishes  with  time.  The  differing  rates 
of  change  in  surface  area  at  different  times  is  a  consequence  primarily  of  particle  collision 
lifetimes  which  change  with  particle  size.  As  time  progresses,  larger  particles  begin  to 
be  comminuted,  supplying  the  smaller  particle  population  at  different  rates  (Sykes  and 
Greenberg  1986). 
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